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ABSTRACT 

The early-type chemically peculiar stars often show strong magnetic fields on their surfaces. These 
magnetic topologies are organized on large scales and are believed to be close to an oblique dipole for 
most of the stars. In a striking exception to this general trend, the helium-strong star HD 37776 shows 
an extraordinary double- wave rotational modulation of the longitudinal magnetic field measurements, 
indicating a topologically complex and, possibly, record strong magnetic field. Here we present a new 
investigation of the magnetic field structure of HD 37776, using both simple geometrical interpretation 
of the longitudinal field curve and detailed modeling of the time-resolved circular polarization line 
profiles with the help of magnetic Doppler imaging technique. We derive a model of the magnetic 
field structure of HD 37776, which reconciles for the first time all magnetic observations available for 
this star. We find that the local surface field strength does not exceed ~ 30 kG, while the overall field 
topology of HD 37776 is dominated by a non-axisymmetric component and represents by far the most 
complex magnetic field configuration found among early-type stars. 

Subject headings: polarization - stars: atmospheres - stars: chemically peculiar - stars: magnetic 
fields - stars: individual: HD 37776 



1. INTRODUCTION 

A sizable fraction of the main sequence A and B stars 
are observed to be strongly magnetic, with typical sur- 
face field str engths between a few hund red G and more 
than 30 kG (jDonati fe Landstreetir2009ft . Unlike the ac- 
tive late-type stars, in which tangled magnetic topolo- 
gies are generated by a contemporary dynamo involving 
differential rotation and turbulence, the early-type stars 
most likely acquire their magnetic fields in the process of 
star formation and the n retain them during the main se- 
quence evolution (e.g.,|Mestcl 2003). To maintain stabil- 
ity, these fossil fields must have a mixed poloidal-toroidal 
geom etry in the stellar interior (Braithwa ite fe NordluncO 
2006), but at the surface the fields are expected to be 
essentially poloidal, with the dipolar contribution domi- 
nating the field geometry. 

Observations of the mean line of sight magnetic com- 
ponent, the so-called longitudinal magnetic field, (B z ), 
generally reveal a s mooth, single- wave variation with 
the stellar rotation (jMathvsl 119911: iBvchkov etaLl 120051: 
lAuriere et al.l T2007) . consistent with a simple, dipolar- 
like structure of the magnetic field topologies in the 
early-type stars. Recent comprehensive investigations of 
the linear and circular polarization spectra suggested the 
presence of small-scale magnetic field structures in some 
chemically peculiar A-type stars (jKochukhov et al.ll2004t 
iKochukhov fe Wadi |2010[ ). At the same time, these 
studies still found a dipolar-like large-scale field organi- 
zation, just as one would expect from simple (B z ) curves 
of these stars. 

Only a few early-type stars show longitudinal 
magnetic field variation deviating significantly from 

a single- wave sinusoid. These stars, HP 37776 

(jThompson fe Landstreetlll985ft . HD 32633 (|Leone et al.l 
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(|Kochukhovl 12006ft . and r Sco (jDonati et al.1 12006ft . 
attracted considerable attention because even a rela- 
tively small distortion of the (B z ) phase curve indicates 
a dominant higher-order multipolar field component. 
This challenges theoretical models to find a satisfactory 
explanation for the existence and apparent stability 
of these rare global stellar magnetic field topologies. 
In this context, investigations of unusual non-dipolar 
fields represent an important step towards achieving a 
more complete theoretical understanding of the hot-star 
magnetism. 

Among the stars with conspicuously non-dipolar lon- 
gitudinal magnetic field variation, HD 37776 (HIP 26742, 
V901 Ori) has the most bizzare (B z ) curve. This young 
He-stro ng B2IV star belongs t o the Orion OB lb asso- 
ciation (Landstrcc t et all I2007T ) and shows a prominent 
rotational vari ation of the brig htness in different photo- 
metric bands (lAdelmarJ [l997ft, spectral line intensities 



(jPedersen fc Thomsenl 119771: iPedersenl 11979ft. and flux 
distribution (Adcl man fc Pvoerl 11985ft due to the com- 
bined effect of a strong magnetic field and an inhomoge- 
neous distribution of chemical elements over the stellar 
surface. Uniquely for an early-type star, its longitudi- 
nal magnetic field curve shows four well-defined extrema 
during each rotation cycle. 

A consistent rotational period of P ro t ~ 1.5387 d fol- 
lows from the variation of the lo ngitudinal magnetic field 
(jThompson fc Land strcct 1985j), the strength of the neu- 
tral He lines dPedersen 19791) . and mu l ti-colo r photom- 
etry (jAdelmanl 11997ft ■ iMikulasek et al.l (|2008ft presented 
a combined analysis of different time-series observations 
obtained for HD 37776 over 31 years. They reported a 
remarkably rapid rotational spin down by 17.7 s during 
this period. 

The first quantitative model aiming to explain the puz- 
zling longitu dinal magnetic fi eld curve of HD 37776 was 
proposed bv iBohlenderl (jl994ft . The (B z ) measurements 
were reproduced with an axisymmetric combination of 
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the dipolar, quadrupolar and octupolar magnetic com- 
ponents with the strengths 3.4, —59, and 44 kG re- 
spectively. The resulting mean surface magnetic field 
strength exceeds the 34 kG fiel d intensity measured for 
Babcock's star ([Babcockl |l960T ). making HD 37776 the 
most magnetized object among all non-degenerate stars. 
However, this conclusion is sensitive to the adopted pa- 
rameterization of the field geometry and assumed incli- 
nation of the stellar rotational axis. The very large pre- 
dicted field strength cannot be confirmed by direct de- 
tection of the Zeeman splitting of spectral lines due to 
t he rapid rotation o f HP 37776. 

iKhokhlova et al.l (|2000D derived magnetic field geome- 
try and surface maps of chemical elements for HD 37776 
with the help of the Doppler- Zeeman imaging technique. 
The magnetic field structure was approximated by a su- 
perposition of non-aligned dipole and quadrupole com- 
ponents. The study employed a simplified treatment of 
the line formation in a magnetic field and was based 
on the circular polarization spectra with a sparse phase 
coverage. Despite these shortcomings, the theoretical 
uvby light curve s calculated for the sur face abundance 
distributions of Khok hlova et al.l (120 00) fit the obser- 
vations very well ferticka et al.l I2007D . On the other 
hand, no definite magnetic field model could be found in 
IKhokhlova et al.l study, with lines of different elements 
giving B d = 10-80 kG, B q = 40-160 kG, and different 
orientations of the magnetic axes. None of these mod- 
els succeeded in reproducing the observed longitudinal 
magnetic field curve. 

The inconclusive results of the previous magnetic field 
studies of HD 37776 and the lack of a definite magnetic 
field geometry model compatible with different types of 
magnetic observations warrants a new investigation of 
this unique star. In this paper, we present a new analy- 
sis of the magnetic field topology of HD 37776 using all 
magnetic field observations available for this star. In Sec- 
tion^ we describe the longitudinal magnetic field mea- 
surements and time-resolved spectropolarimetric obser- 
vations employed in our analysis. Sections [3] and 2] pro- 
vide a description of the multipolar and Doppler imaging 
modeling of the magnetic observations. In Section [5] we 
study the magnetospheric structure of HD 37776 using 
potential field extrapolation. The paper concludes with 
the discussion of our results in Section [6] 

2. OBSERVATIONAL DATA 
2.1. Longitudinal magnetic field measurements 

A stro ng longitudinal magn e tic fie ld in HD 37776 was 
found by iBorra fc Landstreetl (|l!)7!)i. Their seven (B z ) 
measurements showed the presence of a kG-strength vari- 
able magnetic field, which exhibits a peculiar double- 
wave modulation with rotation phase. Longitudinal mag- 
netic field measurements of HP 37776 were continued by 
iThompson fc L andstrcct (1985). All these observations 
were carried out using a photopolarimetric Zeeman an- 
alyzer, tuned to the wings of the H/3 line. In total, 44 
photopolarimetric (B z ) measurements were collected in 
1977-1978 and in 1982-1984, spanning the range from 
-2.2 to +2.5 kG with a typical uncertainty of 300-400 G. 

In 1993, five additional spectropolarimetric longitudi- 
nal field mea s urem ents were obtained for HD 37776 by 
iDonati et al.l (|1997l ). These authors derived the line of 



sight magnetic field component from a mean Stokes V 
profile of 17 metal absorption lines. The errors of these 
measurements are 170-360 G. Despite being of a rather 
different nature compared to the H/3 photopolarimetric 
longitudi nal field diagnostic, the (B z ) data points pub- 
lished by IDonati et al.l ([1997D agree well with the previ- 
ously established longitudinal field phase curve. Taking 
this into account, we combined the photo- and spectropo- 
larimetric (B z ) determinations in the analysis presented 
below. 

We recalculated the rotation phases of the longitudinal 
m agnetic field mea s urem ents using the cubic ephemeris 
of iMikulasek et aT] (|2008l ). which takes into account the 
rapid rotational braking of HD 37776. The updated lon- 
gitudinal magnetic field curve is presented in Figure [TJ 
In this plot the phase 0.0 corresponds to the maximum 
brightness of the star. 

The application of the non-linear ephemeris yields 
a slightly larger scatter of the (B z ) measurements at 
close rotatio n phases compa red to the constant pe- 
riod of, e.g., lAdelman] (|1997l ). Nevertheless, the over- 
all distorted double-wave structure of the longitudinal 
field curve, as well as satisfactory agreement bet ween 
the observations by IThompson fc Landstreetl (|1985f l and 
IDonati et al.l (|1997l) obtained a decade apart, is retained. 

2.2. Circular polarization spectra 

An intensive spectropolarimetric observing campaign 
of HD 37776 took place in 1994-1996 at the 6-m tele- 
scope of the Special Astrophysical Observatory (SAO). 
The star was observed in different wavelength bands with 
the Main Stellar Spectrograph (MSS) equipped with a 
CCD detector and an achromatic circular polarization 
analyzer. Observations of HD 37776 employed the second 
camera of the MSS, giving the effective spectral range of 
5000-7000 A and a dispersion of 14 A/ mm. From Octo- 
ber 1994 to 1996 spectra were recorded using a 530 x 580 
pixel detector, while a 1160 x 1040 CCD was used after 
March 1996. Both detectors provide a spectral resolu- 
tion of R « 10 000, which is adequate for the analysis of 
rapid rotators such as HD 37776. ThAr lamp was used 
as a comparison source. The instrumental profile is close 
to a Gaussian with 2-2.5 pixel sampling of the FWHM. 

The extraction of the MSS Zeeman spectrograms was 
accomplished with the reduction package DECH20 by G. 
Galazutdino\0. This software was also used to perform 
continuum normalization and to determine the wave- 
length solution from the measurement of the centroids 
of ThAr lines in the two orthogonally circularly polar- 
ized comparison spectra. Subsequent steps of the spec- 
tral processing and analysis employ ed the codes writte n 
in the ESO MIDAS environment (|Kudrvavtsevl 12000? ) . 
Further details about spectropolarimetric observations of 
HD 37776 as we ll as reduction and analy sis of these data 
can be found in iRomanvuk et al.l (|1998f ) . 

HD 37776 was observed in three different wavelength 
regions during the spectropolarimetric campaign at SAO. 
By far the most extensive data set of 30 circular polar- 
ization spectra is available for the region containing the 
He I 5876 A line. These spectra have a signal-to-noise ra- 
tio of 300-500 and cover the 5830-5900 A wavelength in- 
terval, which includes the sharp interstellar Na ID lines. 

3 http : //www . gazinur . com/DECH-sof tware .html 
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These features are not expected to show detectable po- 
larization signal, allowing us to verify and fine-tune the 
relative wavelength alignment of the left- and right-hand 
circularly polarized beams in each observation. 

The rotation phases of the SAO spe ctra were calcu- 
lated using the non- linear ephemeris of iMikulasc k et al.1 
(2008), consistently with the longitudinal field data. 
These circular polarization observations provide a dense 
coverage of the rotational cycle of the star. The average 
phase sampling is 0.03 and the largest phase gap is 0.14 
o f the rotation per io d. 

iKhokhlova et al.l (|2000( ) used a subset of 9 spectro- 
grams from the SAO He I 5876 A data in their Doppler- 
Zeeman analysis of HD 37776. This subset provides a 
poor rotational phase coverage, with several phase gaps 
as large as 0.2-0.3. 



Longitudinal magnetic field 
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3. MULTIPOLAR MODELING 

In this section w e aim to revise the multipolar model of 
IBohlenderl ((1994) using the new ephemeris of HD 37776 
and a more realistic value of the incl i nation angle of 
the stellar rotational axis. IBohlenderl <| 1994T ) assumed 
i = 90°, which together with the rotational period found 
in previous studies of HD 37776 and v c sini = 91 kms -1 
determined below (Sect. 01 implies an implausibly small 
radius, R — 2.8-R©, for a main sequence B2 star. Other 
authors reported syste matically larg er radii: R = 3.7Rc?) 
(IShore fc Browdll990h. R, = 4.3i?^ (jZboril et al.l 11997ft . 
R = 3.6i?0 (jLandstreet et al1l2007| ). This range of radii 
corresponds to inclination angles between 40° and 50°. 
Giving more weight to the mo r e rece nt parameter deter- 
mination by lLandstreet et ail (|2007l ). we adopt i = 50°. 
This inclination angle is clo se to i = 45° favored by the 
Doppler imaging analysis of IKhokhlova et all (j2000l ). 

Keeping the inclination angle fixed, we used a least- 
squares optimization algorithm to fit the observed (B z ) 
curve with a low-order, axisymmetric multipolar mag- 
netic field topology. The free parameters of this widely 
used global magnetic field parameterization (jLandstree 
ll988HLandstreet et alJl989HLandstreet fe Mathvsll20Qi ' 
include the angles /3 and 4>, giving the obliquity and az- 
imuth of the magnetic field axis. The field intensity and 
structure are determined by the polar strengths of the 
dipolar quadrupolar B q , and octupolar B a magnetic 
components. 

Our best-fit multipolar model of the field topology 
in HD 37776 is characterized by the reduce chi-square 
= 1.26 and has the following parameters: i = 50°, 
P = 77.1 ± 2.3°, <f> = 184.8 ± 1.7°, B d = -1.0 ± 0.6 kG, 
B q = -4.3 ± 13.6 kG, B = 97.0 ± 7.3 kG. The up- 
per panel in Figure Q] shows the comparison between 
the (B z ) measurements and the model longitudinal field 
curve. Evidently a very strong axisymmetric octupolar 
magnetic field structure is capable of providing a rea- 
sonably good approximation of the complex longitudinal 
magnetic field variation of HD 37776 (Figure [j] upper 
panel). The inferred multipolar field topology is illus- 
trated in more detail in Figure [2j These spherical maps 
demonstrate that the field strength varies between 35 and 
100 kG over the stellar surface. The mean field modulus 
changes from 43 to 49 kG with the stellar rotation phase 
(Figure [1] lower panel). 

Taken at face value, the multipolar field modeling re- 
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Figure 1. Upper panel: comparison of the mean longit udinal 
magne tic field of HP 37776 measur ed by Thompson & Lan dstreeti 
(1985) (circles) and [Donati et aQ l|1997ft (squares) with the fit 
achieved by the axisymmetric multipolar model (dashed curve) and 
by the direct magnetic inversion (solid curve). Lower panel: the 
mean field modulus variation predicted by the multipolar model 
(dashed curve) and by the magnetic DI field map (solid curve). 



suits suggest that HD 37776 is indeed the most magne- 
tized non-degenerate star. 

4. MAGNETIC DOPPLER IMAGING 

We reconstructed the magnetic field topology of 
HD 37776 with the help of the Doppler imaging 
(DI) code Invers IO. This magnetic inversion code 
was described by iPiskunov fe K ochukhov ( 2002|) and 
IKochukhov fe Piskunovl (|2002ft . and applied to the 
analysis of Ap- s tar magneti c field geomet ri es b y 



Kochukhov et al l (120021 I2004D. iFolsom et al 
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InversIO determines the surface distributions of the 
magnetic field vector and chemical abundance from the 
phase-resolved spectropolarimctric observations. The 
code models individual line profiles using a realistic nu- 
merical treatment of the polarized radiative transfer in 
a stellar atmosphere. The surfa ce maps are constrained 
by th e Tikhonov regularization (jPiskunov fc Kochukhovl 
[20021 ). which ensures the uniqueness of the DI reconstruc- 
tion and provides the smoothest solution justified by the 
observational data. 

The absence of the Stokes Q and U data in 
our study might complicate the magnetic field recon- 
struction, leading to some cross-talk between the ra- 
dial and meridional field com ponents at low latitudes 
(|Kochukhov fc Piskunovl [2002f) . Furthermore, as sug- 
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Figure 2. Surface magnetic field structure HD 37776 derived from multipolar fitting of the longitudinal field curve. The star is shown at 
five equidistant rotational phases as indicated above the topmost row of spherical plots. The aspect corresponds to the inclination angle 
i = 50° and vertically oriented rotational axis, a) The distribution of the field strength, with contours of equal magnetic field strength 
plotted every 10 kG. The thick line shows the stellar rotational equator. Rotational axis is indicated by short the vertical bar. b) The 
orientation of the magnetic field vectors. In these vector maps the light arrows show the field vectors pointing outside the stellar surface 
and the dark arrows correspond to the vectors pointing inwards. The arrow length is proportional to the field strength. 



gested by the recent comparison of the magnetic inver- 
sions for the Ap star a 2 CVn using four Stokes param- 
eter observations a nd limiting the input data to circular 
polarization alone (|Kochukhov fc Wa de 2010), the linear 
polarization spectra enable reconstruction of small-scale 
surface magnetic features that cannot be detected with 
Stokes V data. We emphasize that these problems aris- 
ing due to lack of linear polarimetry are not unique to our 
study but apply to all previous Zeeman Doppler imaging 
analyses of cool active st ars and hot massive stars (e.g., 
IDonati et al.ll200l 12006ft . 

For modeling HD 37776 we used a version of our 
magnetic DI code that fits simultaneously the spec- 
tropol arimetric data and the l ongitudinal magnetic field 
curve (jKochukhov et alj|200l . The surface maps were 
sampled on a 1176-element surface grid with an op- 
timal distribution of surface zones (see Figure 5 in 
IPiskunov fc Kochukhov! I2002D . The polarized spec- 
trum synthesis calculations for the He I 5876 A 
line employed an LTE solar metallicity model atmo- 
sphere with T cff = 22000 K and log g = 4.0 from the 
revised Atlas9 gricQ. The effective temperature 
and surface gravity adopted for HD 37776 are consis- 
tent with the stellar para meters determined in pre- 
vious studies of t his star (IGroote &: Kaufmannl 119811: 
Zboril et al]|1997t iHunger fc Groot d 119991: iCi dale et al l 
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Before proceeding to the magnetic inversion, we com- 
pared the predictions of the multipolar magnetic field 
model obtained in the previous section with the observed 
circular polarization profiles of the He I 5876 A line. In 
this calculation the He abundance was adjusted to be 
[He] = \og(N Hc /N R )~\og{N He /N H ) Q = -0.2 to roughly 
match the average strength of this line in the unpolarized 
spectrum. The comparison of observations and multipo- 
lar model predictions is illustrated in Figure[3] Evidently, 
the magnetic field topology inferred from the longitudi- 

4 http: //kurucz .harvard. edu/gr ids .html 



nal field curve alone fails to account for the rotational 
phase dependence of the circular polarization within the 
He I line. In particular, major disagreement between 
the predicted and observed Stokes V profiles is found 
for the rotation phases 0.07-0.29 and 0.61-0.93. Sim- 
ilar significant discrepancies between the details of ob- 
served V profiles of magnetic CP stars and those pre- 
dicted by multipolar m odels have been prev iously found 
for 53 Cam and CrB (iBagnulo et al.l l200l but not for 
78 Vir (|Khalack fc Waddl2006D . " 

In contrast to the unsatisfactory results obtained with 
the multipolar field model, the DI reconstruction of the 
magnetic field geometry yields an excellent fit of the cir- 
cular polarization spectra (Figure[3) and an adequate de- 
scription of the (B z ) curve (Figure [TJ upper panel). The 
resulting magnetic field distribution, constrained by the 
Tikhonov regularization but not limited to a low-order 
axisymmetric multipolar topology, is shown in Figures |4] 
and El The former figure illustrates the spherical surface 
maps of the field strength and orientation while the latter 
one presents the rectangular projections for the radial, 
meridional, and azimuthal components of the magnetic 
field vector. The He abundance distribution obtained 
together with the magnetic maps is shown in Figure 
The He concentration changes from a moderate under- 
abundance with respect to the solar chemical composi- 
tion to an overabundance by up to a factor of 100 inside 
He spots. 

Our magnetic DI reveals a strikingly different picture of 
the magnetic field in HD 37776 compared to the sugges- 
tions of previous multipolar models. The field topology 
is complex and notably non- axisymmetric. There are six 
distinct areas of the positive and negative magnetic po- 
larities, reminiscent of a somewhat distorted £ = 3 spher- 
ical harmonic. At the same time, the DI model yields a 
much weaker surface field strength than commonly as- 
sumed for HD 37776. The field intensity changes from 
» 5 kG to » 30 kG across the stellar surface, which is 
still a lot compared to a typical magnetic B-type star 
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Figure 3. Comparison of the observed Stokes / {left) and V 
(right) line profiles of HD 37776 with the theoretical calculations 
corresponding to the best-fitting magnetic DI map (solid lines) 
and to the axisymmetric multipolar field model shown in Figure [5] 
(dashed lines). The spectra for consecutive rotation phases are 
shifted vertically. The phases are indicated on the right. The bars 
at the bottom give the horizontal (1 A) and vertical (5%) scales 
of the line profiles plots. The horizontal dotted lines identify the 
continuum level for Stokes I and the zero level for Stokes V. 



but certainly less than the field on several most extreme 
members of this group of stars. The mean magnetic 
field modulus of HD 37776 determined from the DI maps 
varies between 13 and 16 kG. This is substantially lower 
than the «46 kG field modulus predicted by the multi- 
polar model (see Figure [U lower panel). 

Aiming to characterize more comprehensively the field 
topology obtained with the magnetic inversion, we ap- 
plied to the final DI maps the spherical harmo nic ex- 
pansion method of lPiskunov fc Kochukhovi l|2002). Both 
poloidal and toroidal spherical harmonic modes with the 
angular degree £= 1-10 and all possible azimuthal num- 
bers m were included in the expansion. The resulting 
expansion coefficients are presented in Figure [6] We 
found that non-axisymmetric £ = 3-4 harmonic compo- 
nents give the largest contribution to the magnetic field 
topology of HD 37776, although a non-negligible power 
is evident for the lower-order multipoles and for the £ 
values up to 5-6. The relative fractions of the poloidal 
and toroidal field components appear to be similar. 

The presence of the significant toroidal field in 
HD 37776 is an interesting and somewhat unexpected 
result. Although the toroidal component frequently 
dominates over th e poloidal one in cool active stars 
(|Donati et al J 12003). it is seldom diagnosed in hot mag- 
netic stars and is never included in the parameterized 
low-order multipolar field models of these stars. Results 
of the present investigation as well as the previous po - 
lari zation line profile stu dies bv lKochukhov et al.l (2004) 
and IDonati et al.1 ()2006D suggest that the toroidal com- 
ponent should not be disregarded in the magnetic field 
analysis of hot stars. 

5. MAGNETOSPHERIC STRUCTURE 

The strong magnetic fields of the early-type stars have 
a prominent effect on their radiatively driven winds and 
circumstellar environment. It is understood that a strong 
enough magnetic field confines the mass loss, produc- 
ing the regions of hot, X-ray emitting shocks and cooler, 
denser clouds w here the plasma is trapped by the mag- 
netic field lines (IBabel fc Montmerle!ll997t IDonati et al.l 
[200lHUd-Doula fc Owockill2002H . The magnetically con- 
fined wind and inhomogeneous magnetospheres mani- 
fest themselves with the c haracteristic var i abilit y of the 
wind-sensitive UV lines (jShore fc Brown! I1990D. emis- 
sion in the hydrogen Balmer li nes (|Groote fc Hunger! 
I1982T) . and photometric indices (jTownsendl 12008( 1 . The 
rotational spin-down, recently de tected in some helium- 
strong stars including HD 37776 (|Mikulasek et al.l 120081 
iTownsend et al.1 [2 01 Of) , is also determined by the inter- 
action between the w ind and the stellar magnetic field 
(|Ud-Doula et al.ll2009h . 

With the exception of the slowly rotating, chemically 
normal early B-type st ar r Sco, which ha s a weak non- 
dipolar magnetic field (|Donati et al.ll2006T) . the previous 
studies of the hot-star magnetospheres were limited to 
the objects with dipolar fields. Our analysis of HD 37776 
probes a previously unexplored parameter range, offer- 
ing an insight into the magnetospheric structure and the 
wind-field interaction for a star possessing a strong mag- 
netic field with a very complex geometry. 

To investigate the extended magnetic field topology of 
HD 37776 we calculated the three-dimensional structure 
of the circumstellar magnetic field using the potential 
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Figure 4. Surface magnetic field structure and He distribution of HD 37776 derived with magnetic Doppler imaging. The star is shown 
at five equidistant rotational phases as indicated above the first and third row of spherical plots. The aspect corresponds to the inclination 
angle i = 50° and vertically oriented rotational axis, a) The distribution of the field strength, with contours of equal magnetic field strength 
plotted every 5 kG. The thick line shows the stellar rotational equator. The rotational axis is indicated by the short vertical bar. b) The 
orientation of the magnetic field vectors. In these vector maps the light arrows show field vectors pointing outward from the stellar surface 
and the dark arrows correspond to the vectors pointing inwards. The arrow length is proportional to the field strength, c) The surface 
distribution of the He abundance, measured relative to the Sun, obtained simultaneously with the magnetic field geometry. The contour 
lines in the spherical abundance maps arc plotted with a step of 0.5 dex. 
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Figure 5. Rectangular projection of the radial, meridional and azimuthal vector components of the magnetic field geometry of HD 37776 
derived in magnetic Doppler imaging. The dotted line shows the lowest visible latitude for the adopted inclination angle i = 50° . The lines 
of equal field strength are shown for the [—30, +30] kG interval, with a 10 kG step. 



field extrapolation technique of Uardine et al.l (|1999f ). In 
this method the magnetospheric field is determined us- 
ing the radial field component of the magnetic DI map 
and setting a spherical source function at the Alfven 
radius, where the field lines are expected to become 
open under the wind pressure. This outer boundary sur- 
face is located where the w ind confinement parameter rj 
(|Ud-Doula fc Owock i 2002), characterizing the ratio be- 
tween the magnetic field energy and the kinetic energy 
of the wind, is on the order of unity. 
Adopting the typical magnetic field strength at 



the stellar surface B+ » 14 . 5 kG , the stellar radius 
JL « AR^ (IShore fc Brownl fl990t IZboril et al.l [HH 
MikulasckltalJ |2008|),~ the average mass-loss rate 
M ~ 10 _9 M Q year " 1 , and the term i nal wind speed 
^oo ~ 10 3 kms" 1 (jKrticka et alj 120061: IMikulasek et al.l 
2008), we find the wind confinement parameter t)~3 x 
10 6 for HD 37776. This yields the Alfven radius R A > 
40.R* for a dipole-dominated surface field and Ra > 10-R* 
for a quadrupolar magnetic topology. 

Despite several uncertainties in our estimate of rj and 
Ra, the results of the potential field extrapolation in 



Magnetic field of the helium-strong star HD 37776 



7 



Poloidal field 



10 
9 
8 
7 

6 - 



6 5 



3 - 
2 - 
1 




4 5 6 7 
/ 

Toroidal field 



10 



10 


i 


i i 


i i i i i i i 


9 








8 








7 








6 








5 








4 








3 








2 










1 















i i i i 


i i i 




1 


2 3 


4 5 6 7 


8 9 10 



/ 

Figure 6. Coefficients for the spherical harmonic expansion of 
the surface magnetic field structure of HD 37776. The top panel 
shows coefficients for the poloidal field, while the bottom panel 
corresponds to the toroidal field component. 



the immediate vicinity of the star are weakly dependent 
on the assumed location of the source surface as long as 
Ra 3> R*- Therefore, we limit the following discussion to 
the inner 2R+ of the stellar magnetosphere. The circum- 
stellar magnetic field topology established in this region 
with the potential field extrapolation procedure is illus- 
trated in Figure [7] In this plot the closed magnetic field 
lines are rendered for four rotation phases. The color 
variation along each field line corresponds to the radial 
field component. 

Not unexpectedly, we find that the magnetospheric 
field structure of HD 37776 is far more complex than a 
dipolar or an axisymmetric low-order multipolar field. 
The extended field topology is dominated by a network 
of the closed magnetic field lines, concentrated along the 
six ridges stretched mostly in the meridional direction. 
The location of the closed field line loops is determined 
by the six distinct areas of the positive and negative field 
polarity on the stellar surface (see Figure [5]) . This non- 
trivial inner magnetospheric field structure should give 



rise to an equally complex rotational phase dependence 
of the parameters tracing the wind inhomogeneities and 
circumstellar matter. 

6. DISCUSSION 

In this paper we presented a new analysis of the mag- 
netic field observations of the unique He-strong star 
HD 37776. The unusually complex magnetic field of this 
object distinguishes it from other early-type stars, which 
have mostly dipolar field topologies. We reanalyzed the 
longitudinal magnetic field measurements of HD 37776, 
finding that a very strong, ~ 100 kG, octupolar field 
is required to reproduce the rotational modulation of 
(B z ) under the assumption of an axisymmetric multi- 
polar topology. However, we found that this multipo- 
lar model is incompatible with the circular polarization 
time-series observations of the He I 5876 A line. Using 
both the longitudinal magnetic field measurements and 
the Stokes V spectra, we derived a new model of the sur- 
face magnetic field topology of HD 37776 with the help 
of magnetic Doppler imaging. This magnetic field struc- 
ture for the first time explained all the different types 
of magnetic observations available for HD 37776. The 
field turns out to be highly complex, bearing no resem- 
blance to a low-order, axisymmetric multipolar geometry 
commonly assumed for early-type stars. T hus, together 
with r Sco, which has a much weaker field ()Donati et al.l 
2006), HD 37776 belongs to an unusual group of the 
early-type stars with conspicuously non-dipolar fields. 

Our investigation of the longitudinal magnetic field 
variations of HD 37776 clearly demonstrated the funda- 
mental limitation of the multipolar field models based 
on fitting the (B z ) curve alone. While the multipolar 
model requires a record-strong field to reproduce the 
(B z ) measurements, magnetic DI yields an order of mag- 
nitude weaker magnetic field intensity, with the aver- 
age field strength w 14.5 kG. This is below the field 
modulus directly mea sured with the Zeeman splitting in 
several Ap/Bp stars dBabcocklll960t iBabel et all 1 19951 
i Mathvs et alJll997t lHubrig et alJl200a iKochukhovl boO 
iRvabchikova et alJl2006t iFrevhammer et al.ll2008| ). 

The He abund ance distribution i n HD 3 7776 was previ- 
ously studied by iKhokhlova et al.l (|2000l ) using a subset 
of data analyzed here and a so mewhat less sophisticate d 
line profile modeling technique (jVasilchenko et aLH l996). 
Their He abundance map is compared with ours in Fig- 
ure [H The He spot distribution derived by Khokhlova 
et al. has a significantly higher contrast and abundance 
range, but the location of the two main overabundance 
features agree reasonably well with our results. On the 
other hand, the fine structure of the He variation across 
the surface of HD 37776 inferred by Khokhlova et al. is 
likely to be spurious. 

Using the magnetic field topology inferred with 
Doppler imaging as a starting point, we applied the po- 
tential field extrapolation method to calculate the ex- 
tended magnetospheric field geometry of HD 37776. The 
field structure above the stellar surface was found to be 
complex and non-axisymmetric, with most of the field 
lines closed within a sphere of 2R+. This magnetic geom- 
etry must correspond to an unusually complicated config- 
uration of the magnetospheric p l asma, which is broadly 
consistent with iShore &: Brown! (|1990T l observations of 
the complex variation of the C IV and Si IV UV wind 
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Figure 7. Magnetospheric structure of HD 37776 inferred with the potential field extrapolation of the surface field determined with 
magnetic DI. The star is shown at four rotation phases and the inclination angle i = 50°. Only the field lines closed within a sphere of 
_R = 2R* are shown here. The color along each field line and across the stellar surface changes according to the value of the radial magnetic 
field component. 



diagnostic lines. 

The existence and stability (at least on the time- 
scale of a few decades covered by the longitudinal field 
measurements) of the very complex magnetic field in 
HD 37776 poses an interesting challenge to the fossil 
field theory, which is usually applied to explain the 
ma gnetism of massive star s. Nu merical simulations 
by iBraithwaite fc Nordlundl (|2006l ) demonstrated that 
both toroidal and poloidal magnetic fields of compara- 
ble strength must be present in the stellar interior to 
maintain the global stability of the fossil field structure. 
On the surface, the field appears to be purely poloidal, 
axisymmetric and dipolar-like. In the context of these re- 
sults, it is difficult to explain complex non-dipolar mag- 
netic fields with a substantial toroidal contribution at the 
stellar surface. 

Subsequent MHD calculations by IBraithwaitel (|2008f ) 
revealed a different class of the stable magnetic equilib- 
rium solutions, consisting of one or several large-scale 
twisted flux tubes lying under the surface of the star. 



Each tube is formed by the toroidal field, with the 
poloidal field lines wrapped around it. The multiple flux 
tube scenario, realized in Braithwaite's simulations if the 
initial unstructured magnetic flux is sufficiently spread 
out inside the star, seems to be applicable to the cases 
of t Sco and HD 37776. However, the theory still has 
difficulties explaining the observed presence of signifi- 
cant surface toroidal field in both stars. In particular, 
it is challenging to understand the origin and nature of 
the current system in the stellar atmosphere required to 

support these toroidal fields. 

The complex equilibrium fields studied by IBraithwaitel 
(2008) evolve on diffusive time-scales, comparable to the 
stellar life-time on the main sequence. The flux tubes 
gradually rise towards the surface and undergo length- 
ening and narrowing due to the loss of toroidal magnetic 
flux. It is generally expected that the diffusive field decay 
proceeds faster for higher-order magnetic components. 
In this case, the magnetic topology of the old stars should 
be closer to dipolar, while the field structure of young 
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Figure 8. Comparison of the He abundance dis tribution derived 
for H D 37776 in our study (top panel) and by Khokhlova ct al. 
lf200GT ) (bottom panel). 



stars , such as HD 37776 (log t = 6.55 vr. lLandstreet et al 
l2007h . could be more complex. However. iBraithwaite 
(2008) presented examples of simpler axisymmetric fields 
evolving into more complex, non-axisymmetric ones. A 
calculation tuned for the specific parameters of HD 37776 
is required to ascertain whether these changes can be ob- 
servable and to what extent the unusually high degree of 
complexity of the field in this star can be related to its 
youth. Nevertheless, we can conclude that the presence 
of a very complex strong magnetic field in a young, early 
B-type star is not inconsistent with the current qualita- 
tive predictions of the fossil field theory. 

We thank Dr. Z. Mikulasek for his help with the ro- 
tational phase calculation for HD 37776. O.K. is a Royal 
Swedish Academy of Sciences Research Fellow supported 
by grants from the Knut and Alice Wallenberg Founda- 
tion and the Swedish Research Council. 
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